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1. INTRODUCTION 

The main advantage of the active suspension system (ASS) is the capability of changing the 
dynamic wheel load (Fdyn) in order to improve the ride comfort for passengers by adjusting the characteristic 
of a suspension system like the stiffness, the damping. The electromagnetic, hydraulic actuators mostly are 
used to generate the active force in the ASS. The control problem for ASS is always more challenging for 
researching and application in the practice because of some reasons [1, 2]: there are many ASS’s states need 
to be controlled while ASS has only one control input, the controlled active force acts contemporaneously on 
both sprung mass and sprung mass, the performance bandwidth of actuator is limited, the limitations of time, 
the response time does not meet the requirement of high velocity, all ASS’s states are not measured by 
sensors, the determining statistic characteristics of road excitation is the practice difficulty. 

There are lots of control methods for ASS including: in the material [3], the authors used the self- 
learning neural networks to adjust controller parameters in order to improve the ride comfort, the complexity 
of neural networks and capabilities of convergence are the disadvantages of this method. Using the Skyhook 
reference model and the adaptive parameter controller based on the Lyapunov method is presented in [4]. 
In [5, 6], authors applied the free-model sliding mode control with time-delay estimation to eliminate 
the nonlinear, uncertainty and noise acting on the system. The modal feedback method used to control ASS is 
written in the material [7], this method control independently in each vibration mode of ASS, the controller 
parameters are turned by an optimal control method. Using second order optimal LQG controller for ASS is 


Journal homepage: http://ijra.iaescore.com 


212 o ISSN:2089-4856 


conducted in the [8], in which the ASS model is rewritten in the linear form. In the article [9], the controller 
based on backstepping method is used to control ASS, in this method control law is designed to back step 
from inner loop to outer loop by using the virtual control inputs. The simple decoupling method for ASS is 
implemented by the authors of [10]. Moreover, there are lots of control methods applied to ASS including the 
Event-triggered control [11], Hoo control method using the Linear matrix inequality, the fuzzy logic control 
method [12], etc. 

The common feature of the above methods is the noise caused by road profile considered to be 
the main noise input of the mathematical model, this model is not concern about that noise which is depend 
on the road profile and velocity of the car. In this paper, we use the general model of ASS, in that model 
the parameters related to the noise caused by road profile and velocity of the car are considered. From 
the practice test, we see that the faster velocity of the car is the small amplitude of the noise is and the faster 
frequency is. This feature makes the noise characteristic acting to the ASS change. Based on that model of 
ASS, the control process for ASS is implemented by two control loops. The inner control loop is designed 
with noises without affecting ASS, this control loop drives the performance of ASS in the transient period. 
The model of closed loop system with inner control loop is used as reference model for the designing 
the outer adaptive control loop. The reference model depends on the velocity of the car. The outer adaptive 
control loop uses the parameter adaptive mechanism adjusted by the velocity of the car and errors between 
ASS’s states and states of the reference model. The adaptive force is generated to compensate for the above 
errors in order to drive that error reaching to zero. The simulation results using the practice road profile 
show that this method is efficient to apply to the practice. The remaining part of this paper is structured as 
follows: Part two mentions the modeling of the suspension system. Part three introduces adaptive control 
for the ASS. Part four shows the simulation results. Finally, part five summarizes some conclusions. 


2. MAIN CONTENTS 
2.1. Modelling of the ASS 

The active suspension depicted in Figure 1 system consists of sprung mass m, and unsprung 
mass Mw, the sprung mass, and unsprung mass are connected by the suspension part including damping 
partc, and spring partd,in parallel [13]. In case of active suspension, the active actuator (electrical, 
electromagnetic, hydric devices driven by electrical signal w(t)) is added to this part to change the natural 
characteristics of suspension part in order to get the desired comfort ride. The connection between 
the unsprung mass and the road surface is the wheel. The wheel is described by wheel damping part c,, and 
wheel spring part dw. The states of the suspension system are deflection, velocity, and accelerator of sprung 
mass Z,,Zs,Z;, unsprung mass Zy,ZyZy, respectively. The z, is the road excitation supported as the main 
input noise of the system, z, depends on the road profile and velocity of the car v(t). 

The performance quality of the suspension system is reflected upon the comfort ride, this quality is 
presented by 3 parameters: the accelerations of the car and the wheel Z,, Z,, respectively, the dynamic wheel 
load F4 and the suspension deflection Zz, — Zw. The dynamic wheel load is calculated as (1) as [14]. 


Fa = Cy (Zr _ Zw) + diy r _ Zw) (1) 





MM 


Figure 1. Physical model of a typical suspension system 


Road profile 
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The dynamic wheel load F4 depends on the mass of the car m,, velocity of the car v(t) and road 
profile. The smaller values of Z,;, Z,,, Fg are, the higher comfort rider of the passengers in the car is Figure 1. 
The road excitation always is supposed as the main input noise, but we realize that the z, not only depends on 
the road profile but also the velocity of the car. For the same road profile, as the velocity of the car changes, 
the stochastic characteristics of the z,will be changed. 

Figure 2 shows the varying of the z, in the case velocity of the car changing from 5 km/h to 45 km/h 
for the same road profile. The road profile has the stochastic characteristics: min = —2.319 cm, 
max = 2.348 cm, mean = 0.0131, median = 0.0131; variance = 0.331678, u = 0.0131, o = 0.576. 
From Figure 3 we see that the velocity of the car is higher, the amplitude of the variation Zz, is lower. 
The maximum of z, is 0.07 ym in the case v(t) = 5 km/h, the maximums of z, are 0.023 m and 0.02 m in 
the cases v(t) = 15 km/h and 45 km/h, respectively. The changing of stochastic characteristics z, by the 
velocity of the car, the road profile will affect strongly to the state estimation process. 

The first order model of the road excitation z, as the function of the road profile and velocity of 

the car depicted in Figure 4 can be written by the equation as [15, 16]: 


dz 
rs = —Av(t)z, + Wz (2) 





where y, v are the optional feedback parameter and velocity of the car, respectively, w,, is the input noise 
caused by the road profile. 
Let consider the state vector of the suspension system described in Figure 1 as: 


Z4 = Zs) Z2 = Za; Z3 = Zw» Z4 = Lites Z5 = Zr) Z6 E Ži (3) 
an T 
Z = |Z Z2 Z3 Z4 Z5 Zel 


Suppose u(t) as the active force, if u(t) = 0 system is called the passive suspension system. The model of 
the suspension system is written by the following equation. 


dz 








= A(v(t))z + Bu(t) + w(t) (4) 
where 
0 1 0 0 0 0 0 
—Cs —d, os. ds 0 0 i Wzs 
Ms Ms Mes Mes Ms, Wazs 
, | 0 0 0 1 0 0 | 0 Wew 
ios Cs ds —Cw-€, —dy—d, Cw- dy Av(t) 0 aim 1 wo Wadzw (5) 
My My My My My E My wae 
0 0 0 0 0 1 0 Wazr 
0 0 0 0 —Av(t) 0 0 
By modifying the (5) as 
dz es 
n A(v(t))z + Bujo (t) + d(t)] (6) 
with 
a _ (u(t) _f 0 z 
B=([B 19@)=("9. dO =|wey]-m = 6 7) 
the system has the input w(t) 
u(t 
aO = (4?) Su =f 0,190 (8) 
n 


Matrix A(v(t)) has six eigenvalues located in the left of the imaginary axis which make 03 pairs of 
eigenvalues that are symmetric with respect to the imaginary axis. Two first of them are fixed and have 
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the negative real part, the last pair depends on the velocity of the car. The velocity of the car increasing 
makes the imaginary part of the third pair of eigenvalues move along the imaginary axis far from the root as 
described in the Figure 3. 


r 


The amplitude of z_(m) 





Time(s) 


Figure 2. The varying of the stochastic characteristics 
for the same road profile by the velocity of the car 
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Figure 3. The locations of pole pairs of the ASS 
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Figure 4. The first order model of road excitation Z, as 
the function of the road profile and velocity of the car 


2.2. Adaptive control for the ASS 
The adaptive control for ASS is designed by two steps as follows. 

Step 1: Design the state feedback controller in order to keep closed loop system without noise 
d(t) = 0 having the desired dynamic. The pole placement is the technique that is used to place the closed 
loop poles of the system in predetermined locations. Based on the desired dynamic wheel load and 
the desired dynamic of Z,,,Z, the predetermined locations of the poles are calculated. The results of 
Step 1 is the state feedback controller which has the gain depending on the velocity of the car and 
the road profile. This is the inner control loop. 

Step 2: Design the adaptive controller in the outer control loop to attenuate the noise caused by the road 
profile and velocity of the car. This control loop uses the reference model to be the model of inner loop 
without noise d(t), the dynamic of the reference model is depended on the velocity of the car and road 
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profile. Using state error between the ASS and reference model, the adaptive mechanism calculates 
the active force to compensate for the noise caused by the road profile and velocity of the car. The 
adaptive mechanism gain is varied road condition and velocity of the car. 


2.2.1. Designing the inner control loop 

To design the controller in the inner control loop we suppose that the velocity of the car is zero, 
so w(t) = O or d (t) = 0,,,. Consider that all states of ASS z(t) are measurable. Therefore, the model of 
ASS can be written as (7): 


dz z 

=, = A(v@)z + Boce) (9) 

So the designing is to determine the matrix R(v(t)) to guarantee that V(t) = 7(t) + R(v(t))z makes 
the model of closed-loop system become 


dz 


ae (A(v(¢)) — BR(v(t))) z+ Bny(t) = A(v(t))z + By) (10) 


and the poles of the dynamic system (8) S1, S2, ...,S, are predetermined locations with respect to the desired 
performance of the ASS [17, 18]. So, the (9) needs to be solved. 


det (s1 — A(v(t)) — BR(w(t))) = (s — s1)(s — 52) (s = s4) (11) 


Note that the matrix R(v(t)) is depended on the velocity of the car v(t). One of the methods to solve (9) is 
the modal method in the material [19]. This method uses the regression calculation process. If mis rank of 
matrix B and /,, A», ..., An are eigenvalues of matrix A then the regression calculation process is implemented 
by n/m steps, on each step the matrix R; is calculated to move m poles of A among /,, 43, ..., An tom new 
poles among S4, So, ..., Sn Of matrix A(v(t)) — BR(v(t)). The matrix R(v(t)) then is calculated by summing 
matrices R;(v(t)). 


2.2.2. Designing the outer adaptive control loop 
Back to the ASS in presence of noise w(t) caused by road profile and velocity of the car (6), with 
the inner control loop the ASS model becomes 


TE = R(E) + Bin + dO] (12) 


in which A is a stable matrix having desired poles $j, S>, ..., Sn. We have the Lyapunov (11) 
A(v(t)) P + PA(v(t)) = -Q (13) 


If symmetric matrix Q is positive definite, (11) always exists a unique positive definite symmetric 
matrix P(v(t)) [20]. 
We use the reference model described by (12): 


Aim = K(v(t))2m + B&O (14) 


Designing the outer adaptive control loop is to guarantee the system (10) tracking to the system (12). 
Suppose d(t) = O(t), we have 


— = A(v(t))Zm + B(C(t) + O(t)) (15) 


So the task of noise attenuation adaptive mechanism is to adjust @(t)in order to make 
the error e(t) = z(t) — Z,,(t) tend to zero, in which 


The algorithm of adaptive control for active suspension systems using pole assign... (Chi Nguyen Van) 


276 o ISSN:2089-4856 


E = AVE )em + BOC) (16) 


The diagram of two control loops for ASS is shown in the Figure 5. 
Use the positive-definite Lyapunov candidate function as follows: 


V(e,0) = eTPe + 07H10 (17) 


in which matrix P is root of Lyapunov function (15) and H is the positive-definite symmetric matrix chosen 


arbitrarily. So to keep ae < 0 ore > 0 for all e + 0, the vector 0 needs to be satisfy (16) 
= dé 
|BPe + H-1—|=0 (18) 
dt 


Therefore, the noise attenuation adaptive mechanism is 


dé = 
— = —HB'P(v(t))e (19) 
dt 
Some remarks: 
Function @&® < 0 is negative-definite with respect to e therefore, @ # d. In another word, the adaptive 


dt 
mechanism (19) is the noise attenuation mechanism, it is not an identification function of noise caused by 


road profile and velocity of the car. 

The velocity of tending to zero of error e is depended on the positive —definite matrix Q. 

Larger the norm of matrix H is the faster error e tend to zero does. Larger eigenvalues of matrix H are the 
larger amplitude of @ is, this causes the oscillation in the system 


dít) 


aii Ne pas = = A(v(t))z + Ble (t)+d Œ) - 


dt 
R(v(t)) 


de = | F 
— = -HBT P(v(t))e (><) 





dim _ A(v(D)em + BEE) 





Figure 5. The diagram of two control loops for ASS: 
inner control loop (blue line), outer adaptive control loop (red line) 


2.2.3. The simulation results using the practical road profile data 
The simulation results below are conducted in two cases that are even road profile and stochastic 
practice road profile. The physical parameters of ASS are: 








Nsec —| 


N N 
1 = 0.2,c, = 8399 (Z) cw = 151176 ( J = 6665 ( E m 500 ( 
m m m 


m 
m, = 194.2 (kg), my = 49(kg) 
a. The simulation results for even road profile 


The initial condition of states is [0;0.2;0;0.2], the amplitude of the event road profile is 0.02 m as 
described in Figure 6. The varying of position, the velocity of the sprung mass and unsprung mass are 
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depicted in the Figure 7. In this figure, the passive case, the pole assign case, and the adaptive case are 
compared with. In the case of pole assign, the number of oscillations of position and velocity of sprung mass 
and unsprung mass is decreased but the amplitude of oscillation is increased in comparison to the passive 
case. In the case of adaptive active control, the amplitude and number of oscillations are both decreased. 
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Figure 6. Event road profile 
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Figure 7. Response of position, velocity of the sprung mass and unsprung mass 
in the case of event road profile 


b. The simulation results for stochastic practice road profile 

Probability density function of the stochastic practice road profile (as in Figure 8) is: 
min = —2.319cm, max = 2.348cm, mean = 0.0131, median = 0.0131; variance = 0.331678, u = 
0.0131. The simulation results are conducted for three velocities of the car 5 km/h, 40 km/h and 
80 km/h. Figures 9 and 10 are depicted the simulation oscillation responses of ASS, the dynamic wheel 
load and active force, respectively. In these figures, the active control responses are compared to the passive 
case. Similarly, Figures 11 and 12 show the simulation results in the case that the velocity of the car is 
40 km/h, Figures 13 and 14 plot the results for 80 km/h. 

At 5 km/h of the velocity of the car, we see that the oscillation and amplitude of oscillation 
the sprung mass is effectively decreased as shown in the Figure 9. The active control force acts to the system 
clearly at 0.5 s from the beginning. Efficiency of the oscillation decrease of the car is better than the wheel. 
Reducing car oscillation is more efficient than reducing the wheel oscillation. The dynamic load varies in 
range +300 N, the active control force ranges from -600 N to 400 N. 
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Figure 8. Stochastic practice road profile 
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and unsprung mass for cases: passive and control force at 40 km/h (velocity of the car) 
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adaptive control at 80 km/h (velocity of the car) 


At car speed of 40 km/h with active adaptive control, the car oscillation decrease is quite good 
in comparison to the passive case. For the first 0.5 seconds of the adaptive control process, the dynamic 
wheel load varies in the range -300 N to 200 N, and after that the dynamic wheel load decreases, it ranges 
in +100 N. The active control force varies from 400 N to 400 N for the first 0.5 seconds and it varies in 
the range +200 N for after that. At car speed of 80 km/h, the simulation results are the same at 40 km/h. 
Therefore, by using the cascade adaptive control for the ASS, the oscillations of the car and the wheel are 
decreased effectively, it is the same case with the dynamic load. 


3. CONCLUSION 

The paper introduces the cascade adaptive control method for an ASS. The inner control loop uses 
the pole placement technique which is used to place the closed loop poles of the system without noise in 
predetermined locations in order to change the dynamic characteristics of ASS to desired performance 
requirements. The outer control loop uses the adaptive control strategy with the adaptive mechanism to 
compensate for the oscillation of the car caused by the road profile and velocity of the car acting to ASS. The 
simulation results using even road profile and stochastic practice road profile show that the capability of 
oscillation decrease for the car and the wheel of ASS is quite efficient. The dynamic wheel load is decreased 
so the ride comfort of passengers is better. 


ACKNOWLEDGEMENTS 
This research is supported financially by Thai Nguyen University of Technology, TNUT. 


REFERENCES 

[1] D. Karnopp, “Theoretical Limitations in Active Vehicle Suspensions,” Vehicle System Dynamics, vol. 15, no. 1, 
pp. 41-54, 1986. 

[2] M. C. Smith, “Achievable Dynamic Response for Automotive Active Suspensions,” Vehicle System Dynamics, 
vol. 24, no. 1, pp. 1-33, 1995. 

[3] J. Lin and R. Lian, “Intelligent Control of Active Suspension Systems,” JEEE Transactions on Industrial 
Electronics, vol. 58, no. 2, pp. 618-628, 2011. 


The algorithm of adaptive control for active suspension systems using pole assign... (Chi Nguyen Van) 


[4] 


[9] 


[10] 


g ISSN:2089-4856 


A. Turnip, et al., “Controller design for active suspension system based on skyhook reference model,” in 2015 
International Conference on Technology, Informatics, Management, Engineering & Environment (TIME-E), 2015, 
pp. 147-151. 

V. S. Deshpande, M. Bhaskara, and S. B. Phadke, “Sliding mode control of active suspension systems using 
a disturbance observer,” in 2012 12th International Workshop on Variable Structure Systems, 2012, pp. 70-75. 

W. Jue and Z. Jing, “Model-free tracking control for vehicle active suspension systems,” in 2015 34th Chinese 
Control Conference (CCC), 2015, pp. 8067-8072. 

F. Braghin, F. Resta, and E. Sabbioni, “A modal control for active/semi-active suspension systems,” in 2007 
IEEE/ASME international conference on advanced intelligent mechatronics, 2007, pp. 1-6. 

B. Zhang, M. Fan, and F. Miao, “Optimal Control of Vehicle Active Suspension Systems with Actuator Delay,” in 
2007 IEEE International Conference on Control and Automation, 2007, pp. 2257-2261. 

A. Alleyne and J. K. Hedrick, “Nonlinear Control of a Quarter Car Active Suspension,” in 7992 American Control 
Conference, 1992, pp. 21-25. 

G. Tang, H. Lin, and H. Su, “Decoupling vibration control for active suspension systems,” in 2017 29th Chinese 
Control and Decision Conference (CCDC), 2017, pp. 6504-6509. 

K. Bansal and P. Mukhija, “Event-triggered control of vehicle active suspension systems,” in 2078 Indian Control 
Conference (ICC), 2018, pp. 178-183. 

S.-j. Liu, Z.-h. Huang, and Y.-z. Chen, “Automobile active suspension system with fuzzy control,” Journal of 
Central South University of Technology, vol. 11, no. 2, pp. 206-209, 2004. 

C. N. Van, “State Estimation Based on Sigma Point Kalman Filter for Suspension System in Presence of Road 
Excitation Influenced by Velocity of the Car,” Journal of Control Science and Engineering, vol. 2019, 2019. 

G. Koch and T. Kloiber, “Driving State Adaptive Control of an Active Vehicle Suspension System,” in JEEE 
Transactions on Control Systems Technology, vol. 22, no. 1, pp. 44-57, 2014. 

D. Hrovat, “Survey of advanced suspension developments and related optimal control applications,” Automatica, 
vol. 33, no. 10, pp. 1781-1817, 1997. 

G. Koch, et al., “A nonlinear estimator concept for active vehicle suspension control,” in Proceedings of the 2010 
American Control Conference, 2010, pp. 4576-4581. 

K. Ogata, Modern Control Engineering. Prentice Hall PTR, 2001. 

C. N. Van and P. N. Doan, “Adaptive tracking control based on disturbance attenuation and ISS stabilization of 
Euler-Lagrange nonlinear systems in the presence of uncertainty and input noise,” in 20/1 2nd International 
Conference on Artificial Intelligence, Management Science and Electronic Commerce (AIMSEC), 2011, pp. 3698- 
3701. 

G. Yao, et al., “Modal computing method of output feedback control gain matrix in pole assignment of vibration 
system,” in 2010 International Conference on Mechanic Automation and Control Engineering, 2010, pp. 2243- 
2247. 

G. Kitagawa, “An algorithm for solving the matrix equation X = FXF T + S,” International Journal of Control, vol. 
25, no. 5, pp. 745-753, 1977. 


Int J Rob & Autom, Vol. 9, No. 4, December 2020 : 271 — 280 


